n = 10;
alpha = 1;
sigma2 = [0.82, 0.91, 0.16, 0.18, 0.28, 0.2, 0.04, 0.07, 0.36, 0.97]';
pmin = 0;
pmax = 1.5;
G = [13.0 0.47 0.84 0.95 0.2 0.05 0.8 0.09 0.48 0.71; 0.54 29.0 0.91 0.8 0.87 0.48 0.51 0.95 0.06 0.44; 0.26 0.35 27.0 0.12 0.65 0.67 0.39 0.78 0.47 0.28; 0.91 0.77 0.72 11.0 0.59 0.84 0.48 0.18 0.38 0.34; 0.17 0.04 0.12 0.08 21.0 0.17 0.07 0.08 0.98 0.24; 0.66 0.27 0.77 0.78 0.74 29.0 0.48 0.68 0.22 0.07; 0.58 0.07 0.8 0.1 0.55 0.78 18.0 0.93 0.86 0.11; 0.87 0.16 0.04 0.84 0.48 0.51 0.45 11.0 0.07 0.89; 0.97 0.61 0.48 0.18 0.6 0.62 0.64 0.65 29.0 0.51; 0.77 0.14 0.9 0.31 0.79 0.31 0.67 0.6 0.85 23.0];
